Dimensional cross-over and charge order in half-doped Manganites and Cobaltites 
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We propose a generic model for understanding the effect of quenched disorder on charge ordering 
in half-doped Manganese- and Cobalt-Oxides with different crystal structures. Current experimental 
observations (Table are explained in the light of the global phase diagram of the model (Fig. 0. 



In the past decade much theoretical work was devoted 
to elucidating the nature and origin of charge ordering 
(CO) in the doped transition metal oxides, particularly 
in manganiteSjWhich are of great experimental and prac- 
tical interest 0, 0, El • Much attention was focused on 
studying the interplay of charge correlations with orbital 
and spin degrees of freedom which are often important 
for the CO Surprisingly, however, no coherent ap- 
proach to understanding the effect of quenched disorder, 
such as introduced by a random distribution of dopant 
ions, on charge order (distinguished from orbital order) 
has been proposed so far. 

Dopant disorder is a source of an unavoidable random 
electrostatic potential which couples linearly to charge 
density fluctuations 0] . In a strongly correlated electron 
system close to criticality such perturbation is of cru- 
cial importance. We propose that essential account for 
this random potential can be formulated in terms of a 
generic model Ijljl. mapping the CO interaction on an ef- 
fective Ising problem. We construct the phase diagram of 
this model, Fig. and find that it adequately describes 
the essential features of the CO structure in half-doped 
manganese and cobalt oxides. In particular, a dimen- 
sional cross-over in our model gives a natural explana- 
tion to the puzzling experimental observation that CO 
in the layered perovskites is short-range 0, 0J3Li while it 
is long-range in the pseudo-cubic materials 0,0, El ■ 

At half-doping, the ground state of Co oxides and of 
many isostructural Mn oxides manifests a checkerboard- 
type planar charge modulation accompanied with a Jahn- 
Teller (JT) distortion of the M0 6 (M = Mn, Co, etc.) 
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FIG. 1: Charge order and JT distortion at half doping: ab 
plane is tiled with cells {r} containing two cobalt (manganese) 
ions, (a)-(d). Crystal structure of pseudo-cubic perovskite 
material Pro.sSro.sMnOs (e) and layered Lai.sSro.5Co04 (f). 



octahedrashowninFig.Hl HE 0,0, SEES- Theelec - 
trons are highly localized in a charge-ordered state with 
alternating Co 2+ and Co 3+ valence (Mn 4+ and Mn 3+ in 
manganites), and the in-plane 2_ ions move towards 
the higher- valence M ions. To be explicit, we present 
the model construction in terms of a Lai.sSro.5Co04 0] 
which has no superimposed further symmetry breaking 
due to orbital order. 

As a general framework, we argue that hopping dy- 
namics, orbital, magnetic, local Coulomb and Jahn- Teller 
interactions can all be integrated out into an effective 
charge-ordering interaction which then competes with 
the random charge potential introduced by the dopant 
ions. A perfectly ordered state has a two- fold degener- 
acy associated with the choice of sublattice occupied by 
the higher/lower valence ions. We define a local Ising 
variable, r (r) = ±1, to denote the two alternative ionic 
configurations, Fig. [He). A single-domain CO state cor- 
responds to a ferromagnetic order of r (r) . Topological 
disorder is introduced by domain walls as shown in Fig. 
^d). In principle, the model can be enriched by intro- 
ducing the 'vacancies', r (r) = 0, which represent cells 
with even number of electrons, but we shall exclude such 
possibilities in the analysis and discussion of this paper. 

Because of the Jahn- Teller distortion, there is an elas- 
tic strain energy cost for nearest neighbor cells with op- 
posite r values, Fig. ^d). This strain energy can be de- 
scribed by an effective short-range ferromagnetic interac- 
tion J(r- t')t (r)r(r') Without a loss of general- 
ity we may approximate J (r — r') by a nearest-neighbor 
coupling. The potential introduced by dopant ions ran- 
domly favors one or the other configuration of a two- 
ion cell, and therefore effectively acts as a random field 
h (r) r (r) . Importantly, depending on the crystal struc- 
ture the effective coupling may be anisotropic, J c ^ J a b- 
Hence, for the charge order, the effect of elastic crystal 
strain and random dopant distribution is captured by an 
effective anisotropic random- ficld-Ising- model (RFIM), 
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where (rr') , and (rr') denote nearest-neighbor cells in 
the a6-plane and along the c-axis respectively, see Fig. Q] 
(c), (f). To our knowledge, all previous analytical and nu- 
merical studies focused on the isotropic 3D RFIM. Here 
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FIG. 2: Charge order phase diagram from the anisotropic ran- 
dom field Ising model Q . J c / J a (, and J a b / A parameterize the 
anisotropy of the effective coupling and the relative strength 
of the disorder potential from dopant ions, respectively. 

we make the first attempt aimed at a qualitative and 
quantitative understanding of the anisotropic RFIM QJ. 

Random field alters both lower critical dimension 0] 
and critical exponents |14| of an Ising system. In partic- 
ular, for weak random field, A/J-^i 1, the isotropic 3D 
RFIM orders at a finite temperature T c > 0, while the 
2D RFIM has no long range order (LRO) even at T = 
(J and A = {rms of h r } parameterize the Ising coupling 
and disorder potential, respectively). This prompts an 
interesting question: what is the fate of a 3D anisotropic 
random field Ising model (ARFIM) of QJ? The simplest 
way of connecting the renormalization group (RG) flows 
leads us to the phase diagram shown in Fig. [2 

To connect our model analysis with the experiment, 
we summarize the charge order correlation lengths in the 
afe-plane, £ a ;,, and along the c-axis, £ c , measured in sev- 
eral representative perovskite oxides in Tabic ^ Cor- 
relations in the pseudo-cubic systems I-III (Fig. ^ (e)) 
arc resolution-limited (typically this means £ > 2000 A 
00); an d CO is apparently long-range. On the other 
hand, in the layered systems IV-V (Fig. Q (f)) the CO 
is finite-range, with highly anisotropic 3D correlations 
£ab ^ £o which implies J a b S> J c - Large in-plane cor- 
relation length £ a b 3> 1 indicates weak planar disorder, 
A <C J a b- Therefore, these materials probe the most 
interesting and complicated region of the phase diagram. 

Because the isotropic 31? RFIM is long-range ordered 

TABLE I: Correlation lengths in terms of number of unit cells 
in each direction (i.e. £ c in units of c ~ 12.5 A, and £ a b in 
units of ay/2 « 5.44 A, as implied by Fig. 0. 
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for A < A 3 C D w 2.27 J [ig, it is clear that for A < J c < 
J a b the system is in the regime of weak 3D disorder and 
has a 31? LRO. Thus, we examine an intriguing case, 



J c < A < J ab 



(2) 



where the c-axis coupling is in the regime of strong dis- 
order, while the a6-planes are in the weak disorder limit. 

For J c — the system decouples into independent 2D 
RFIM planes, and the ground state is disordered. This 
form of disordered ground state is perturbatively sta- 
ble. For an infinitesimal inter-plane coupling J c /A <C 1, 
the afr-planes remain uncorrelated on all length scales 
(c-axis correlation length £ c < 1), and each plane may 
still be treated as an effectively independent 2D RFIM. 
Therefore, a critical line of phase transitions exists in 
the (J a b/ A, J c / Jab) phase diagram, separating the disor- 
dered and 3D-LRO phases (solid line in Fig. 2). With in- 
creasing J c / Jab on the disordered side we enter a critical 
regime, where the ground state is still disordered but with 
significant and highly anisotropic correlations in all spa- 
cial directions. In the isotropic 2D and 3D RFIM the 
ground state correlation length scales with k = A/ J as 



i 2D « exp(i K - 2 ), and £ 



3D 



(ft — k c ) (for ft > k c ), 



respectively. What are the correlation lengths in a highly 
anisotropic 3D RFIM JTJ), with J c « A < J ab l 

The phase diagram, Fig. instructs that upon a RG 
transformation model JD scales either to a 3D LRO fixed 
point, or towards a 3D strong disorder fixed point. Cur- 
rently, there is no generally satisfying analytical deriva- 
tion of the RG equations for the RFIM for d < 4. How- 
ever, because d = 2 is the lower critical dimension of the 
RFIM, 3D random field critical fixed point should be ac- 
cessible perturbatively in 2 +e expansion, Therefore, 
we derive perturbative RG equations for the anisotropic 
3D RFIM starting from a semi-phenomenological 
scaling "ansatz" for the 2D RFIM |I5ll7|. These equa- 
tions, which naturally account for the RG flow of the 
coupling anisotropy, allow us to derive and analyze all 
qualitative features of the phase diagram of Fig. 

For the isotropic <i-dimensional RFIM a simple scal- 
ing analysis of Imry-Ma, j J (I) - Jl^ 1 ; A (I) - VFA j, 
(TsL . predicts the relative disorder k = y to scale as 
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The marginality of k (I) in 2D is associated with d = 2 
being the lower critical dimension. It is broken by sub- 
dominant contributions from the interface roughening at 
the domain boundaries which lead to a disordered ground 
state For weak disorder, ft <C 1, the correlation 

length is £ = Acxp | + ift^ 2 }, with a w 1 (clearly, this 
expression fails for strong disorder, where £ < 1). By 
demanding that ft (I) scales consistently with the cor- 
relation length, £ (I) — £/Z, and with the expression 
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£ (/) = Aexp { + -^k 2 (/)}, we obtain a scaling equation 

K ~ 2 (l) = K ~ 2 [1 - K 2 a In (I)) , (4) 

i.e. scaling of the relative disorder in 2D is governed 
by the logarithmic corrections to Eq. The disor- 

dered ground state is manifested by the flow of k^ 1 (I) = 
— > 0. In the context of the model QJ, equation Q 
describes the RG flow for J c = 0. 

To proceed with deriving the qualitative features of the 
RG flows shown in Fig. [3 for < J c < J a b we note that 
the anisotropy of interactions is, in fact, an irrelevant 
perturbation at the 3D critical point. This is reflected 
by the fact that all stable fixed points are located on 
the isotropic line in the phase diagram, at J c /J a b = 1- 
At the mean-field level (and without random fields) the 
anisotropy of the interactions can be removed by appro- 
priate anisotropic re-scaling of the coordinates [ljj . Here 
we implement this idea by devising an anisotropic real- 
space RG transformation, effectively coarse-graining to 
anisotropic blocks l a b x l a b x l c . We derive the RG equa- 
tions along the lines of 2 + s expansion, and demonstrate 
explicitly the flow towards the isotropic fixed points. 

We parametrize the anisotropy as a — , a; > 0, 

and define the RG transformation through an infinitesi- 
mal anisotropic re-scaling, {d(lnl a b) = d(ln/) ; d(ln/ c ) = 
a(l)d(\nl)}. This determines the shape of the coarse- 
grained blocks to be gradually varying, flowing from the 
anisotropic towards isotropic. We use the Imry-Ma ap- 
proximation 01 to derive the lowest order contributions 
to the flow equations for the model parameters, and sup- 
plement it by the appropriate higher order correction to 
reproduce the scaling equation (@J for the 2D RFIM case, 
a = 0. The resulting perturbative RG equations are 
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We obtain an unstable 2D fixed point at {a* — 0, k* 2 = 
0}, and stable isotropic 3D fixed points at a* = 1. The 
3D critical point is at {a* = 1, k* 2 = KZ = a* /a}. With 
k c m 2.27 from numerical simulations 15], the extrapola- 
tion of our perturbative equations to a = 1 would imply 
a ~ 0.2, which is too small. This indicates that while 
qualitatively correct, and valid in the perturbative re- 
gion of parameters, RG equations ©,© are inadequate 
for quantitative estimates. 

It is straightforward to integrate Eq. JBJ) and obtain 
ail)' 1 = 1 + (a(l)- 1 - l)r x , the scaling flow of the 
anisotropy. In principle, the phase transition line can 
be determined by a further numerical integration. The 
important asymptotes can be derived analytically, 

o.(k 2 ) » Co e~~^ , 

a(K 2 )«l + Ci (k 2 -k 2 ), 1 
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FIG. 3: Block transformation procedure 

and show that neglecting O (a 2 ) term in Eq. JSJ is indeed 
justified for a < 1, but not for a ~ 1. 

To obtain quantitative estimates of the anisotropic cor- 
relation lengths as a function of J Cl J a b and A, we use 
an alternative approach. In the absence of the exact 
RG equations for the 3D ARFIM, we resort to a real- 
space rescaling procedure where we employ the Migdal- 
Kadanoff bond- moving technique |18| combined with the 
scaling ansatz lf4fl. We adopt conventional approxima- 
tions in ignoring the generation of longer range interac- 
tions and using the "majority rule" |l8j for block spin 
variable. Hence, we follow only the transformation of ef- 
fective nearest neighbor block interactions J a b(l), Jc(J) 
and effective disorder parameter A(Z). 

Our strategy is to map the anisotropic RFIM onto 
a "solvable" isotropic 3D RFIM, whose known proper- 
ties enable us to derive the estimates for £ a b (J a b, Jc, A) 
and £ c (J a b, Jo A). Because the correlation length scales 
as £ (I) = a bare anisotropic 3D model with finite 
S,ab > £,c > 1 is mapped onto an effective isotropic 3D 
model with £ a f, = £ c when the coarse-grained block size 
Lb x l ab x l c is chosen with l ab = £ a b/Cc and l c = 1. 
In other words, the cross-over from quasi-2D to isotropic 
3D scaling is intuitively implemented in two steps. First, 
a 2D transformation to blocks of size l a b x l„j x 1, and 
then a usual isotropic 3D scaling (see Fig. 3). 

Performing the Migdal-Kadanoff transformation to 
lab x l a b x 1 blocks we first move the J c bonds ^| and 
thus create an effective block interaction J c (l a b) = Jclab' 
Fig. 3B. Now that the c-axis coupling is removed, each 
planar l a b x l a b section can be integrated out in a 2D 
fashion using J3J. This approximation is equivalent to 
assuming that fluctuations on length scale l a b < £ab/£c 
in afo-planes are uncorrelated in the c-axis direction. In 
addition, we keep the Imry-Ma approximation for scal- 
ing of the disorder A (l a b) under 2D block transformation 
|17j . Indeed, a perturbative expression for the surface 
tension per unit length on the scale I in 2D RFIM [lfij ]. 
E(Z) = J [l — (tk 2 ln(/)l , indicates that logarithmic cor- 
rections to K (I) come from the effective weakening of the 
inter-block interaction J (I) as domain boundaries adjust 
to the random field. We thus arrive at the final form of 
our approximate scaling equations for block l a b X l a b X 1 , 
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For l a b = £,ab/£,c, the correlation lengths in the trans- 
formed model are isotropic £ a b (lab = Cab/Cc) = £c- Since 
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the disorder remains isotropic, it must be that the 
rescaled interactions are also isotropic. Thus, by im- 
posing a self-consistency condition, J a b(lab = £ab/£c) = 
Jc (£ab/£c)> we obtain a mapping of the bare anisotropic 
RFIM onto isotropic 3D RFIM with effective inter- 
actions 



Jab (U = W&) = Jc (W£c) - [f Jc , (12) 

and A(£ afc /£ c ) = (jr) A - Usin g GD^dD.GS) and 

expression for the correlation length of the 3D RFIM, 
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Substituting £ 
in the phase diagram 



1, we obtain the crossover line £ c ^ 
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With k c « 2.27, A « 0.1 and f=s 1.37 from numerical 
simulations, [l5|. equations (|13|) . l|14l) allow us to retrodict 
the effective model parameters. For Lao.5Sri.5Co04 (Ta- 
ble 1, V) we find A « 0.85, and J-f- > 53, indeed satisfy- 
ing the relations J c <C A < J a (, under which our approx- 
imations are valid. For Lao.sSri.sMnC^ (Table 1, IV) 
-A- pa 0.64 and ^ « 49. In both cases the anisotropy of 
the bare coupling parameters is remarkably high! Yet, we 
think it is quite realistic for the layered systems. Because 
neighbor afe-planes are shifted so that M06 octahedra in 
one plane fit in-between those in the other (see Fig. 1(f)), 
and the inter-plane spacing is rather large, breathing- 
type distortions of the octahedra accompanying the CO 
are weakly coupled between the planes. On the other 
hand, octahedra within each plane form a corner-sharing 
network, so the coupling of in-plane distortions is strong. 
In the pseudocubic perovskite structure, the M06 octa- 
hedra share apical oxygens along the c-axis as well (see 
Fig. 1(e)), and crystal fields are more isotropic. Note 
that, for a system to have a LRO groundstate, interac- 
tions need not be isotropic. In fact, the phase diagram 
instructs us that only very strongly anisotropic systems 
will be disordered if the disorder potential is weak (as ex- 
emplified by Eq. J7J), and the above estimates of J a b/J c 
in the layered compounds). 

In conclusion, we argued on very general grounds that 
charge order observed in half-doped manganites and their 
isostructural relatives is described by the anisotropic 
(J c 7^ Jab) 3D random field Ising model We con- 
structed the schematic phase diagram of this model, Fig. 
2, and supported its intuitive structure by perturbative 
RG equations JSJl , © . We also obtained the quantitative 



estimates of the effective model parameters as a func- 
tion of the measured correlation lengths using a Migdal- 
Kadanoff block transformation scheme combined with a 
phcnomenological scaling equation (0} for a 2D RFIM. 

Long-range charge order is observed in the pseudo- 
cubic materials at low temperatures. In contrast, a disor- 
dered ground state with anisotropic correlations is found 
in the layered compounds. This disparity is naturally ex- 
plained in our model: because of the strong anisotropy of 
the layered materials these systems reside in the different 
regions of the phase diagram. The anisotrpic correlation 
lengths measured in Lao.sSri.5Co04 and Lao.5Sri.5Mn04 
corroborate our conclusion. We note that a disordered 
3D RFIM has well known features which may be veri- 
fied experimentally. Furthermore, our analysis may be 
applied to thin films, where even in pseudocubic com- 
pounds correlations in one direction are limited by the 
film thickness. The latter, together with the model pa- 
rameters, will then determine the CO finite correlation 
length within the film. Finally, we note that charge dis- 
order in the form of dislocations, as in Fig. Efd), has 
inevitable effects on the superimposed orbital and spin 
orders which have thus far not been elaborated. 
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